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Molecular dynamics simulations were conducted to investigate the dynamic properties of melts of nonconcate- 
nated ring polymers and compared to melts of linear polymers. The longest rings were composed of iV = 1600 
monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. 
The ring melts were found to diffuse faster than their linear counterparts, with both architectures approxi- 
mately obeying a. D ^ N^"^-^ scaling law for large N . The mean-square displacement of the center-of-mass 
of the rings follows a sub-diffusive behavior for times and distances beyond the ring extension {B?g), neither 
compatible with the Rouse nor the reptation model. The rings relax stress much faster than linear polymers 
and the zero-shear viscosity was found to vary as ryo ~ jv^ '**'^'^ which is much weaker than the N^''^ behavior 
of linear chains, not matching any commonly known model for polymer dynamics when compared to the 
observed mean-square displacements. These findings are discussed in view of the conformational properties 
of the rings presented in the preceding papei^S. [DOI: 10.1063/1.3587138] 



I. INTRODUCTION 

Understanding the static and dynamic properties of 
ring polymer melts is one of the remaining challenges in 
polymer science. Unlike linear polymers, the topolog- 
ical constraints for ring or cyclic polymers are perma- 
nent and this affects both their static and dynamic prop- 
erties. For linear polymers the topological constraints 
imposed by the non-crossability of the chains (entangle- 
ments) force each chain to diffuse along its own coarse- 
grained backbone, the so-called primitive path, and this 
is well described by the reptation model of de Gennes 
and Edwards. For branched systems strands have to fold 
back in order to find a new conformation without cross- 
ing any other chain, resulting in an exponential growth 
of the longest relaxation time due to the entropy bar- 
rier of ©(strand length) between different stated^^^^. A 
number of simulation and experimental results confirm 
this conceptP^. For both linear and branched polymer 
melts, it is the free chain ends that make the known re- 
laxation mechanisms possible. However, in the case of 
rings there are no free ends. This raises a number of 
unanswered questions regarding the motion and stress 
relaxation of ring polymer melts. 

In the companion papeJ^' we discussed the difficulties 
of understanding the equilibrium (static) properties of a 
melt of nonconcatenated ring polymers. We pointed out 
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that this system presents a formidable challenge to the 
current polymer theories. It is of considerable current 
experimental intereslP^t^o]^ and it is believed to have a 
significant projection on the problem of chromatin struc- 
ture and dynamics in chromosomes. We refer the reader 
to the detailed discussion of these issues in the Introduc- 
tion section of our papei'i^ the large bibliography of rele- 
vant references is also in that paper. The difficulty of the 
melt of rings problem is such that even the equilibrium 
structure of nonconcatenated rings is highly non-trivial. 
We have computed a large number of equilibrium char- 
acteristics of this system in the previous papei'i'. 

In the present paper, the second in a series of two, 
we investigate the dynamic properties of ring polymer 
melts using molecular dynamics (MD) simulation. To 
our knowledge the only previous MD studies of ring 
polymer melts are by Hur et al.E2 and Tsolou et al.'22l 
who considered polyethylene for short chain lengths (i.e., 
^ ^ 3A^e, linear), and Hur et al.l^for chains up to about 
llAi'e, linear, wherc iVe, linear IS the entanglement chain 
length of linear polymers at the same density. Here we 
consider rings and linear melts composed of much longer 
chains. The number of monomers per ring ranges from 
100 to 1600 or 37Ve, linear to 57iVe, linear • Measuring chain 
lengths in terms of A^e, linear provides a basis to compare 
rather different simulational as well as experimental stud- 
ies. The linear chain melts were simulated using the same 
model as the rings with chain lengths varying from 100 
to 800 monomers per chain. Such highly entangled melts 
have very long relaxation times and a highly parallelized 
simulation code had to be used. For the longest rings as 
many as 2048 IBM Blue Gene/P cores were used for a 
single simulation. 
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All details about the model and the setup of the sys- 
tems can be found in the preceding paper, where the 
very same simulations have been analyzed in order to 
understand the static properties of nonconcatenated ring 
polymers. All properties reported here are exclusively 
taken from the fully equilibrated part of the data. The 
simulation model is briefly presented in the next section. 
This is followed by a section on the dynamics as well as 
a primitive path analysi^^H of the ring systems. A com- 
prehensive discussion of the results is given in Section [IV| 
with the key findings of the work and future challenges 
presented in the Conclusions section. 

II. SIMULATION METHODOLOGY 

Based on the Kremer-Grest (KG) modeP^l, all beads 
interact via a shifted Lennard- Jones potential (WGA 
or Weeks-Ghandler- Anderson potential) with a cutoff of 

= 2^/^ cr. Nearest neighbor beads along the chain in 
addition interact via the finitely extensible nonlinear elas- 
tic (FENE) potential, and chain stiffness is introduced 
by an angular potentiaL^i^The model parameters are the 
same for all ring and linear polymer simulations reported 
in this work. 

The natural time scale of the WGA potential is r = 
cr-y/m/e, where m is the mass of a monomer. For the 
simulations we use T = 1.0 e/ks and p — 0.85 a^^. This 
model has been frequently used for simulation studies of 
polymer melts and solutions, so that one can refer to am- 
ple information throughout the discussion of the present, 
new results. Also, the linear chain data presented in 
this work is all new. The ring melts were simulated us- 
ing 200 polymers of length 100 to 1600 monomers per 
chain. The linear systems were composed of 2500 chains 
of length N = 100, 250 chains of TV = 200, and 400 
chains of length N = 400 and 800. An important pa- 
rameter for our discussion is the entanglement length for 
linear polymers in a melt. For the model employed here 
this parameter was determinec^ to be 

A^e.lincar = 28 ± 1 . (1) 

For more details of the model, the setup and the nota- 
tion and software used see the preceding papei'^. 

III. DYNAMICS RESULTS 

A. Mean-Square Displacements of Monomers and Chains 

The dynamics of ring polymers pose special problems 
for a theoretical description. Based on the non-trivial 
structural properties found in the preceding paper we 
have to anticipate dynamic properties not observed in 
other polymer systems. While for short rings one can ex- 
pect that the Rouse model might describe the dynamics 



and stress relaxation fairly well, this problem for longer 
rings is completely open. From the static properties one 
sees that for large N the rings exhibit an TV-independent 
intrachain density, Pscif , resulting in a scaling of the over- 
all extension as a globule, while at the same time there 
is still significant interpenetration between rings. For a 
dense, unknotted and nonconcatenated globule one could 
expect a Rouse-like dynamics for the beads on short time 
scales and an overall diffusion coupled to amoeba-like 
shape fluc tuatio ns, similar to segregated chains in two 
dimensions i^^ l ^^ l However, the mutual interpenetration of 
rings results in a stronger interchain coupling of the mo- 
tion. Whether there are entanglement effects as in linear 
polymers is not clear since there are no free ends to rep- 
tate. Kapnistos et al.ESl and Milner and NewhalP^ have 
developed models for stress relaxation based on a lattice 
animal picture. This however requires entropically un- 
favorable double-stranded conformations, which are only 
somewhat close to the conformations corresponding to 
the lattice animal scaling regime of Fig. 2 in the previ- 
ous papei'^. Thus since the expected behavior is quite 
unclear, we present here a rather comprehensive first 
analysis of single and collective ring polymer dynamic 
properties. The motion of the rings and linear systems is 
characterized by the various mean-square displacements 
(MSD), and for the rings these results do not provide 
a consistent picture in terms of the known theoretical 
models. 

The first quantities to look at are the different mean- 
square displacements of monomers and whole chains as 
a function of time. The MSD averaged over all the 
monomers of a chain is gi(t), the MSD of the monomers 
with respect to the center-of-mass of the chain is 52(^)7 
and the MSD of the center-of-mass, rcM, of the chains 
is 53 (i): 

5i(<) = (|r.W-r,(0)|^), (2) 

52(<) ^{\r^it)~ remit) - r,(0) + rcM(0)|') , (3) 

g3{t) = {\rcM{t)-rcM{Q)\'). (4) 

In Eqs. (2) and (3), is the position of monomer number 
i in space. 

For linear polymers gi{t) and 32 (0 ^-^^ averaged over 
the middle five beads to reduce chain end effectJ^^. In 
all cases the diffusion of the center-of-mass of the whole 
simulation system is subtracted. Fig. [T] shows gi{t) and 
g3(t)N for the ring and linear systems. Rouse and repta- 
tion theory provide a reasonab ly go od explanation for the 
behavior of the linear system a^ l ^" ! For the mean-square 
displacement of the monomers, giit), the reptation con- 
cept predicts in its asymptotic form** a sequence of power 
laws of gi{t) - t^/'^,t^/^,t^/'^,t^. For short times the 
standard Rouse behavior dominates. After reaching the 
Rouse relaxation time Tg of an entanglement length Nf, 
the constraints imposed by the surrounding chains cause 
the beads to follow a Rouse motion along the random 
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FIG. 1. (a) gi{t) and (b) g-i(t)N versus time for the linear systems. These quantities are shown for the rings in (c) and (d). 
The expected scaling behavior from Rouse and reptation theory is shown along with the t^^'^ regime for the rings. 



walk contour of the reptation tube, resulting in a t^/^ 
power law. After the Rouse relaxation time tr ~ N"^ the 
whole chain diffuses along the tube contour, resulting in 
the second t^/^ regime followed by the free diffusion in 
space after the chain leaves the tube, in the ideal model at 
time ^ -^"^/-/Ve, linear- This intermediate slowing down 
of ^ ^^^"^ comes along with a sub-diffusive regime 
g^it) ~ i^/^ for the displacement of the chain center- 
of-mass. It should however be noted that the expected 
(73 (i) ~ t in the early Rouse regime typically reduces to 
a power law due to correlation hole effects. 

Fig. [T^,b gives the results for the linear polymer 
melts for chains up to iV = 800. With an estimated 
-^e, linear = 28 ± 1 wc cxpcct a wcll developed reptation- 
like behavior for the longer chains. We observe a Rouse- 
like behavior for the shorter chains and a clear interme- 
diate regime in gi{t) with an exponent close to t^'^^ sep- 
arating the two t^/'^ regimes in gi{t) for the long chains. 
There is the well known very smooth crossover, making 
it very difficult to precisely identify the crossover times 
for the different regimes. Fig. [Tjo shows g^{t)N versus 
time for the linear systems. As expected from the Rouse 
description there is an immediate crossover towards free 
diffusion for the shorter chains whereas the longer ones 
exhibit a well-defined t^/"^ regime before crossing over into 
free diffusion (with some small deviations for short times 
for N — 800). In all cases the free diffusion of the center- 
of-mass of the chains sets in long before g^ (t) reaches the 
radius of gyration squared, i.e., at about /2.5 for 

N = 200 and (Rl) /3.0 for N = 800. Though there are 
still many details to be understood in the context of the 



dynamics of linear polymers, the Rouse and reptation 
theories give a very good semi-quantitative description. 

While for linear polymer melts the basic dynamic con- 
cepts are fairly clear, the situation for ring melts is strik- 
ingly different. Fig. [ij: shows gi{t) for the rings. All 
five data sets follow the same curve during early times. 
Then, independent of N we observe a significant slowing 
down in the mean-square displacement from the initial 
i^/^ scaling to a power law very close to the t^/^ ex- 
pected for perfect reptation of linear polymers. Note 
however that the rings cannot entangle in the classi- 
cal sense and thus the rings are not expected to show 
reptation-like behavior. A scaling regime of gi{t) ~ t" '^^ 
has been reported for melts of polyethylene ringp2] of 
length 11 A^e. linear which agrees with the present findings. 
After the approximate t^/* regime we do not observe a 
second t ' regime in gi (t) at a value of 0{N^/^) (much 
smaller than the radius of gyration squared) , but rather 
a direct crossover toward the diffusion of the whole chain. 
From Fig. [l]: we observe this crossover at a displacement 

significantly larger than (i?^)^^^, something not observed 
in any other polymeric systeml The crossover times are 
around t w 7.6 x 10* t, 2.4 x 10^ t, 1.5 x 10^ r, 8.2 x 10^ r 
and gi{t) w 150 a^, 200 a^, 330 ct^, 450 cr^ for N = 100, 
200, 400 and 800, respectively. For the N = 1600 case, 
no onset of a crossover towards diffusion is seen for up 
to 2.0 X 10'^ r. In Fig. [iji, g3{t)N is shown for the rings 
along with a reference line with a slope of 3/4. This scal- 
ing behavior is seen to provide reasonable agreement with 
the center-of-mass motion of all five ring systems and for 
the N = 400 system it holds for several decades. The 



4 



10'' 



(a) N = 200 




10* 



10-^ 



b 102 



10^ 



10" 



: (b) N = 800 


1 ' ' 


1 ' ' 


1 ' E 


























t : 



]3 10* 


10^ 


10" 


10^ 



t/r 



FIG. 2. The mean-square displacements are shown as a func- 
tion of time for the rings and linear systems with A*' = 200 
and 800. The inner monomers are used for the calculation of 
gi{t) and g2{t) for the linear chains. The scaling laws from 
Rouse and reptation theory are also shown. 
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power law reminds one of a correlation hole effect 
for Rouse chains, rather than any reptation phenomenon. 
The crossover towards free diffusion is in good agreement 
to what is observed from gi{t). However, the crossover to 
the free diffusion time scaling of occurs at values sig- 
nificantly larger than the chain extensions (cf. Table 2 
of the preceding paper), similar to what is seen for gi{t). 
This is in contrast to the linear systems which show a 
scaling after moving even somewhat less than their own 
size. For the longest rings the fully diffusive regime is 
not reached within the time frame of the simulation. For 
N = 1600 eventually we observe an even smaller effec- 
tive power law, deviating from the t^^^ towards a mini- 
mal power around The mean-square displacement 
data for the rings and linear systems is available in the 
supplementary materiaP^. 

Another striking observation from Fig. [TJi is that 
g3{t)N for short times does not become independent of 
N. It displays a small but still clearly visible ring length 
dependence, indicating a delayed onset of the constraints 
imposed by surrounding rings with increasing ring size. 



Such phenomena have not been observed in a polymeric 
system so far and suggests long-range ring-ring correla- 
tions, even though the chains cannot be entangled in the 
classical sense. The slight increase of g3{t)N for short 
times and increasing N is indicative of the fact that crum- 
pled globules, which is the presumed state of nonconcate- 
nated rings in the melt, are characterized by a rather deep 
correlation hol^^l]. 

To gain additional insight into these correlations we ex- 
amined the dynamics of the exchange of neighbors. The 
neighbors of a given chain are defined as the other chains 

1 /2 

whose center-of-mass is within a distance {R^) of the 
center-of-mass of the reference chain. In the preceding 
pa.Tpe^, we denoted the average equilibrium value of this 

1/2 

quantity as Ki{r) with r — (i?e) the root-mean- 

square spanning distance between beads N/2 apart along 
the ring. At any given time each chain has a specific set 
of neighbors with the average number of chains in the 
set being Ki, which is 11.6 for N = 100 and 16.7 for 
N = 1600. We computed the number of remaining (un- 
exchanged) neighbors after the reference ring moved its 
own root- mean-square gyration radius. This number was 
found to increase with increasing chain length from 5.5 
for N = 100 to 7.3 for N = 1600. We also computed how 
far the reference chain moved before all neighbors but one 
were exchanged. The calculation was performed using a 
large number of time origins and the root-mean-square 
displacement was found to increase with increasing chain 
length from about 2.5 (i?^^^/^ for N = 100 to signifi- 

1/2 

cantly more than 3(i?^) for N = 800. It is interesting 
to note that these displacements are roughly the same as 
those required for the systems to reach the fully diffusive 
regime. 

Even though the ring polymers display a remarkable 
and so far unseen slowing down compared to a Rouse- 
like motion, they diffuse significantly faster than their 
linear counterparts. This is illustrated in Fig. [2] where 
the MSD for two chain lengths are compared directly for 
both melts of rings and of linear chains. As expected 
.92 (i) and 53 (i) of the rings cross around 2 (^R^y For 
the longer rings, however, 33 (i) at this crossing point has 
not yet reached the fully diffusive regime. For the lin- 
ear polymers the entanglement effects are much stronger, 
and the crossing of g2{t) and 53 (i) is hardly reached. It 
should however be noted that this crossover for the linear 
chains is below 2 (i?^) and occurs in a regime where 33 (i) 
already displays the standard diffusion behavior. The 
self-diffusion coefficients, as determined by the long-time 
behavior of 33 (i), are shown in Fig. [sjfor the rings and 
linear systems. The rings diffuse faster than the linear 
systems for all chain lengths considered. This is in agree- 
ment with previous simulation resultJ^^MH. And while 
early experimental paperJi^EIl report that the two ar- 
chitectures have similar diffusivities, recent experiments 
cast doubt on these findingji^l. Reptation theory predicts 
for linear polymers D ~ for N > iVejinear; however, 
in practice the observed scalin^^H jg ]y-2.3±o.i ^^^^ ^^lis 
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FIG. 3. Self-diffusion coefficient multiplied by A'^ versus chain 
length for the rings and linear melts. The line through the ring 
data was fit to the values for the N — 400 and 800 systems 
while for the linear data all points were used. The arrow 
beneath the point for the largest ring system indicates that 
its value is expected to decrease somewhat until the system 
reaches the fully diffusive regime. The range of this shift 
(length of arrow) is based on the extrapolation of g3{t)/t to 
infinite time. With the exception of the A'^ = 1600 ring case, 
all error bars are smaller than the size of the symbols. 
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FIG. 4. (a) Linear chain with N = 800 shown every 12000 r 
for 60 conformations or 0.28rij_iincar and (b) ring with N = 
1600 shown every 10000 r for 60 conformations or 0.23t_b, rings. 
Chain snapshots at later times are drawn darker in color with 
the final chain conformation in black. 



is found here to within the error bars. The data for the 
ring systems are almost parallel (in log-log scale), but 
shifted to significantly faster diffusion. Using the values 
from the N = 400 and 800 cases one finds D ~ N^^-^ for 
the rings. This is a stronger scaling than reported pre- 
viously for shorter rings. One reason for this is that the 
present study uses longer chains with iV/iVejincar ~ 57 
for N = 1600. Despite a simulation time of more than 
3.9 X 10^ T, the longest ring system did not reach the 
fully diffusive regime. The arrow in Fig. [3] indicates the 
extrapolated final value which will be slightly lower than 
the data point. Because the linear systems follow the 
reptation picture for chains as short as iV = 100, the ef- 
fective entanglement length for the rings must be higher 
than that of the linear chains, if the data were to be inter- 
preted in terms of some as yet unknown generalization of 
the reptation concept. The corresponding entanglement 
length for the rings, A'e, rings, will be estimated in Sec- 
tion IIIB This is all to be compared to a visualization 
of the chain motion, which significantly differ between 
the two architectures. Fig. |4] compares the motion of 
a linear chain to that of a ring for times less than th. 
The linear chain is found to be largely confined to a tube 
with the end monomers having more freedom than the 
inner monomers. This is in sharp contrast to the ring in 
Fig. [4]d which is found to trace out a roughly spherical 
volume. 

The slow crossover towards free diffusion of the rings 
is quite striking. For melts of strictly two-dimensional 
polymers it is known that the chains segregate and form 
2-d globular structures. The diffusion is mostly due to 



amoeba-like shape fiuctuations, which do not require the 
complete reorganization of the internal winding of the 
chain. In a similar way we can check whether the diffu- 
sion time of a ring, the correlation time of the radius of 
gyration squared and the internal contour diffusion follow 
the same time dependency or not. If a and b are span- 
ning vectors whose tails are separated by N/A monomers 
(see inset of Fig. 5) then c = a x b. In practice, the 
points forming the head and tails of the spanning vectors 
are averaged over five adjacent monomers. The correla- 
tion time of the autocorrelation function Ccc{t) provides 
a time scale for the internal structural rearrangements of 
the rings. Correlation times are computed by numerically 
integrating the normalized Ccc(i) from zero to infinity. 
These times are compared to the diffusion times and the 
correlation time of the radius of gyration squared in Fig. 
[5] The longer rings are found to undergo substantial in- 
ternal rearrangements much faster than they diffuse their 
own size or relax the overall extension. Not only is the 
amplitude different, but also the exponent in the power 
law dependence of N on Trciax is significantly lower and 
approaches a value close to 2.2. Thus the overall diffu- 
sion decouples from the local internal rearrangement, as 
implied by the decay of Ccc{t)- This strongly suggests 
that the slowing down in the mean-square displacement 
of the overall chain, 53 (t), is due to ring-ring coupling, 
which then also dominates the overall relaxation of R^. 
For N > 400 the correlation times of are found to be 
significantly larger than the internal reorganization times 
and they closely follow the time scaling of the overall dif- 
fusion of the chains. This implies that internal structural 
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FIG. 5. The diffusive and conformational relaxation times 
divided by as a function of chain length. Inset: The a 
vector is formed between monomers i and i + N/2 while b 
is formed from i + N/A and i + ZN/A. With c = a x 6, 
the time scale for single-chain structural relaxation is found 
by integrating the normalized autocorrelation function of c. 
Despite a very long simulation time it was not possible to 
accurately determine the correlation time of Cjj2fl2 it) for the 
N = 1600 ring system. The error bars for the relaxation 
times of the ring systems with < 1600 are smaller than the 
symbol sizes. The line through the linear chain data is a fit 
while the Treia^/N ~ N line is for reference only. 



rearrangement takes place with a fairly small change in 
the overall extension of the chain. 

If the internal reorganization time and the diffusion 
time are that different, this might also have consequences 
for the viscosity of the melts. The shear relaxation mod- 
ulus is 

GMt) = j^{aMt)a^p{0)), (5) 

where a'Q^(t) is the pre-averaged stresP^. In Eq. [5] 
a and /3 are Cartesian indices with a ^ (3, and V is 
the volume of the simulation box. G{t) is the average 
of Gai3{t) with a/3 = xy, xz, yz. The shear relaxation 
modulus is shown in Fig. [6^,b for the linear and ring 
systems. The linear systems show the characteristic be- 
havior. That is, a Rouse-like decay at early times fol- 
lowed by a plateau region and finally exponential decay. 
For the rings the behavior is quite different. The ring 
systems show G{t) ~ followed by exponential de- 

cay. Specifically, the exponent for the stress relaxation 
tends to increase from —0.5 to —0.45 with increasing N . 
The functional form and exponent for G{t) are nearly the 
same as those suggested by Kapnistos et The ring 
melts are found to relax stress much more rapidly than 
the linear systems. There is no sign of a plateau even for 
the largest rings, suggesting that a standard tube model 
certainly is not applicable. The difference in noise be- 
tween the two data sets arises from differences in averag- 
ing schemes and system sizes (cf . Table 1 of the preceding 
paper). For the linear chains, the stress was computed 




FIG. 6. Stress relaxation for the (a) linear and (b) ring sys- 
tems. The Rouse scaling is shown for reference in (b). G{t) 
was computed as described in the main text but additional av- 
eraging was performed to make the curves more discernable. 
The error bars for G{t) increase with increasing t. 



every other step for 25 steps and averaged for each 100 
step interval. For the rings the stress was computed at 
every step and averaged over each 1000 step interval or 
10 r. 

Results for the zero-shear viscosity obtained from the 
Green-Kubo relatiorP^, 770 = G{t)dt, are shown in 
Fig. [7] The linear systems show a N^-"^ scaling which 
is in agreement with previous simulatiorPsl and exper- 
imental results.'^ Even though the error bars become 
quite large for the longer rings and for TV = 1600 we can 
only give a rough estimate of 770, the rings show a much 
weaker dependence on chain length with r]o ~ jyi-4±o.2^ 
For comparison, the Rouse model predicts a linear A'^- 
dependence. This is another result which requires a 
deeper understanding. It appears that stress relaxation, 
which governs the viscosity, follows a time behavior much 
closer to the internal reorganization of the rings rather 
than the overall diffusion and relaxation of R^. The vis- 
cosity measurements for the rings and linear systems have 
been confirmed by nonequilibrium molecular dynamics 
simulation^Sll These results will appear in an upcoming 
papei!^. The zero-shear viscosities are given for the rings 
and linear systems in Table |l] 

As an additional check, a simulation was performed 
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FIG. 7. Zero-shear viscosity of the rings and Hnear melts 
versus A'^. The error bars for the A'^ = 400 ring case were 
obtained from four independent simulations while the other 
points were obtained from single simulations. While it is diffi- 
cult to estimate the error bars for the TV — 1600 case those of 
the 800 case are expected to be smaller than twice the symbol 
size. Note the arrow above the data point for the A'^ = 1600 
system. 
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TABLE I. Transport coefficients for the rings and linear sys- 
tems. For the same value of TV the rings are found to have a 
higher diffusivity and a lower viscosity. D and rjo are given in 
units of a^/r and ra"^ /e, respectively. The diffusion coefficient 
and the viscosity for the TV — 1600 ring case are upper and 
lower bounds, respectively, since both have not settled down 
completely within our analysis window. This is partially due 
to the extensive averaging needed. The error bars for rjo for 
the TV = 400 ring case were obtained from four independent 
simulations. 



on the TV = 800 ring system using a dissipative par- 
ticle dynamic^^SEIl (DPD) thermostalP^ instead of the 
Langevin thermostat. The same friction coefficient was 
used and Tc was taken as 2i/6 cr. The simulation was 
run for 8 x 10^ t and it showed the same behavior for 
the mean-square displacements indicating that hydrody- 
namics or momentum conservation are unimportant for 
these systems on this time scale. 



B. Primitive Path Analysis 

Even though ring polymers in a melt cannot entan- 
gle in the classical sense, the slowing down of the over- 
all chain diffusion as well as the properties of gi (t) and 



the slow exchange of neighbors suggest that neighboring 
rings must be somehow coupled in order to prevent fast 
free diffusion. Thus an adapted primitive path analysis 
(PPA) was conducted to elucidate the nature and extent 
of possible quasi-entanglements in the ring systems. The 
procedure for linear chain^^^l -^^^s modified for this pur- 
pose. Starting with an equilibrated configuration, a pair 
of opposite monomers {i and i + TV/2) of each ring poly- 
mer are chosen at random to be fixed in space throughout 
the PPA. With the angular and nonbonded intrachain in- 
teractions switched off, the temperature was lowered to 
0.05 e/fcs and the dynamics were carried out until equi- 
librium was reached. Without fixed monomers the rings 
collapse to points. At least three cases were run for each 
value of TV. 

Fig. [8]shows the last configuration from each PPA sim- 
ulation for TV = 100, 400 and 1600. For the TV = 100 case, 
many of the chains have collapsed to double-folded lines 
with an average contour length, Lpp, of 2(i?e). Rings 
that have not completely collapsed are mutually con- 
strained by one or more neighboring chains. As TV in- 
creases, the equilibrium fraction of double-folded lines 
decreases dramatically. For instance, for the TV = 1600 
case shown in Fig. [Sj; all the chains are mutually quasi- 
entangled with not a single double-folded line. It can also 
be seen for increasing chain length that there is a greater 
tendency for rings to thread or pass through one another. 
While the chains were initially segregated at the start of 
the simulation, at equilibrium they are found to form a 
highly connected melt, even though the volume of each 
ring is only shared by a small number of other rings. 

While somewhat questionable here, the standard argu- 
ment to deduce an entanglement length for linear chains 
may be applied to the rings. For the linear chains the 

entanglement length is (^-Rg unoar) /^pp^ which is the 
number of monomers per Kuhn segment of the primi- 
tive path. For the rings this expression is modified by 
recognizing that a ring should be treated as two lin- 
ear chains each of length TV/2 with common end points. 
In Fig. |9] we plot TVg^iings against chain length where 
iVe,ri„gs = (TV/2)(i?2_^i„g,)/(ipp/2)2. For the linear 
chains one finds TVejinear ~ 28, which is independent of 
TV for sufficiently long chains. The entanglement num- 
ber for the rings is found to vary with chain length with 
the longest rings giving approximately 77. This value, 
however, still seems to slowly grow, which is not totally 
surprising in view of the early time values of g3{t)N as 
shown in Fig. [ijl and the slow crossover to an asymptotic 
value of Ki (see preceding papei^. The average contour 
length for each ring system is given in the inset to Fig. 
|9] The final configuration from the PPA for the TV = 800 
linear system is shown in Fig. 10 for comparison. The 
entanglement density of the linear system with TV = 800 
is higher than that of the TV = 1600 ring system. Our 
computed value of TVe^uncar for this system agrees with 
the previously reported value of 28 ± 1-^-. 

These results are to be compared with the fact, found 
in the previous paper, that the pressure increase due to 
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(c) Rings N = 1600 



(a) Rings N = 100 (b) Rings N = 400 

FIG. 8. Final configurations from a primitive patli analysis of the ring systems: (a) A'^ = 100, (b) 400 and (c) 1600. 




FIG. 9. N^. rings as a function of chain length. For our model 
A*'e, linear = 28± 1. The iuset gives the average contour length 
versus A'^. 



Linear N = 800 




FIG. 10. Same as Fig. [8] except for a linear system with 
M = 400 and iV = 800. The linear chains with iV = 800 are 
far more entangled than the A*' — 1600 rings (cf. Fig. [sf. 



the entropy loss associated with the rings squeezing each 
other in the melt was too small to be observed in our sim- 
ulations. If the concept of entanglement length iVe^rings 
makes sense for rings then it should be possible to esti- 
mate this pressure increase as l/iVe ^ings times the refer- 
ence pressure. Given the above estimates of iVe. rings we 
do not expect to be able to observe such a small difference 
in pressures. 



IV. DISCUSSION 

In the previous papei^S we investigated the structural 
properties of a melt of nonconcatenated rings. It was 
found that the overall conformations have a number of 
peculiar properties and are by no means trivial. These 
peculiar structural properties of the long ring polymer 
melts have significant consequences for the dynamics. 
The diffusion of linear chains is explained by the Rouse 
and reptation theories as demonstrated in the present 
work. These theories make use of the fact that linear 
chains have free ends but ring polymers do not have free 
ends so the application of these theories, and especially 
reptation theory, is questionable. However, the Rouse 
model can be solved for a Gaussian ring^^i. The diffusion 
times from Rouse theory for the linear and ring mod- 
els are given in Table [ll] Even for the smallest value of 
TV only fair agreement is found between the simulation 
results and those predicted by the theory for both ar- 
chitectures. The disagreement grows with increasing N 
as the systems become more and more entangled, a phe- 
nomena which is completely ignored by the theory. The 
Rouse theory prediction for the viscosity of the rings is 
in better but also not perfect agreement with the scaling 
law obtained from simulation. 

The peculiar conformational properties of the rings are 
expected to lead to new dynamic properties. The ring 
melt simulations give a number of unexpected and not 
fully understood results. For linear polymers the overall 
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Rings 


Linear 




Rouse Theory 


Simulation 


Rouse Theory 


Simulation 


N 


td/10" giiro) 


td/IO" giiro) 




td/10" giiro) 


100 


0.16 


62.1 


0.07 


47 


0.32 94.7 


0.48 78 


200 


0.65 


125.1 


0.39 


91 


1.29 178.3 


6.77 182 


400 


2.59 


249.8 


2.38 


164 


5.12 355.0 


59.7 283 


800 


10.4 


500.5 


20.0 


292 


20.7 714.1 


559.4 


1600 


41.4 


998.6 


> 105.5 




82.9 1429.2 





TABLE II. Comparison of the diffusion times and mean-square displacement of the monomers between Rouse theorj'^'^ and 
the simulation results. The diffusion times are computed using td ~ {R1)/6D. For the Rouse model D — kT/NC, for both 
architectures while {^) is Nl^l/G for the linear and Nl^l/Vi for the rings. To compute td from the simulation data, D is 
obtained from Table ifand {Rg) from Table II of the previous paper^. The MD model is mapped to the Rouse model with 
Ikl ~ 2.71 (T^ and ^ = 43 r^^. For the linear chain Rouse model (/i(td) was computed using the middle of the chain (i.e., 
n — N/2). Td and giro) are given in units of r and , respectively. 



stress relaxation and the crossover towards free diffusion 
of the beads is linked to a motion distance of the order 
of the radius of gyration, even though the center-of-niass 
motion, 33 (i), follows a power law at much shorter 
displacements in comparison to the rings. This is due 
to the fact that any relaxation is connected to the poly- 
mer diffusing out of its original tube. Such a tube con- 
cept certainly cannot be applied for rings, at least not 
directly. Also the strong connection of diffusion distance 
to stress relaxation does not apply any more here with- 
out significant alterations. In short, the dynamics can be 
summarized by the following observations: 

• sub-diffusive behavior of 173 (i) ~ up to distances 
2-3 times (^g) with y « 0.75, and a diffusion time 
To ^ N^, X > 3 (for N ~ 1600 we even find a 
slowing down 53 (t) ^ , y 0.65, for larger times), 

• the "^1/4" regime in gi{t) also extending 2-3 times 

• diffusion coefficient D{N) ^ N""^-^ depends on N 
in a similar way as in entangled linear polymers, 
however the prefactor is approximately 7 times 
larger for the rings with N — 400 and 800 than 
the corresponding linear polymers, 

• zero-shear viscosity scaling 770 ^ ^i-4±o.2^ which 
is much more suggestive of Rouse-like behavior 
TV) than reptation (~ TV^-^), is in accord with an 
internal reorganization time ©(iV^-^^"-^). 

The first three points seem to contradict the fourth 
which suggests that the overall motion of the chains and 
the stress relaxation, at least to within our accuracy, 
are somewhat decoupled. In the following we propose 
a mechanism which can account for this. Also, the fur- 
ther slowing down of gait) for N — 1600 is required in 
order to prevent an early, unphysical crossing of gi (t) and 
gait) (see below). 

At short times (< < t^) we expect the motion of the 
beads, 51 (i), to follow the standard Rouse t^^'^ behavior, 
and as shown in Fig. [l] this is indeed found. This is 



in contrast to 53 (i) which does not follow the standard 
Rouse behavior as discussed above. Once the entangle- 
ment time is reached, a slowing down at a motion ampli- 
tude of the order of the tube diameter ~ ^^-^g linear 
is expected. This is observed for both the linear and 
ring polymers. Then one expects for linear polymers the 
Rouse motion along the tube, leading to the well-known 
^1/4 regime and then, after the Rouse time, the diffusion 
out of the tube, resulting in the second t^^'^ regime, until 
after the disentanglement time Td ^ N^'"^ where the free 
diffusion dominates. For rings we expect a similar initial 
slowing down due to the motion along the short branches 
of the lattice animal-like structure (shorter rings) or of 
the crumpled globule (large rings). It appears, however, 
that this slowed down motion and the diffusion along 
the chain contour itself, as suggested by the decay of the 
autocorrelation function of the cross product of the span- 
ning vectors, cf. Fig. [5] create a unique power law very 
close to t^/"^ over a very long time scale and distances 
exceeding the gyration radius of the rings. Then there is 
a crossover towards free diffusion without an additional 
clearly-defined power law in between, which is in contrast 
to the behavior of linear chains. 

Although we do not have a consistent theoretical ex- 
planation for the observed dynamics, we can formulate an 
argument relating some of our observations to the others. 
Specifically, the slowed down sub-diffusive motion of the 
beads requires some reduction in the effective exponent 
of the center-of-mass motion as well, otherwise gi it) and 
gait) would meet or even cross (meaning that the coil 
as a whole would move much further than a monomer, 
which is obviously impossible) at a time even well be- 
low the Rouse time. Let us assume 53 (t) ~ t"^ for the 
"chain diffusion" at intermediate times and 51 (i) ~ t"^ 
for one monomer diffusion. We might write for a time 
dependent "diffusion constant" £)(gi(t))p2 = g-^^it). As- 
suming that Digiit)) ^ giit)^^^^^, namely proportional 
to the number of beads carried along at a distance g 
within an object of fractal dimension df — 3, we arrive 
at the relation = (5/2)ail^. This relation appears to 
be consistent with our simulation data. Indeed, careful 
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inspection of the data for gi (t) indicates that the small- 
est slope is close to 0.28, systematically decaying with 
increasing chain length to slightly larger than 1/4. The 
smallest slope in 53 (i) for N = 1600 is close to 0.7, which 
is in good agreement with the above relation between ai 
and 0:3. 

Also the amplitude of 33 (t), where the crossover to- 
wards free diffusion occurs needs some more attention. 
For instance, as shown in Fig. fl| for the N — 400 ring 
case the transition to 53 (t) ~ t ai long times occurs at 
roughly 7.1 x 10^ r where gsit) = 170 ct^ or 3.2(7?^). For 
the linear case at the transition one finds g3{t) = 170 
or 0.9(i?g). For the ring systems with N = 400 and 
800 (and extrapolated for TV = 1600) one then finds 
D ^ N^'^-'^. All previous simulation papers have re- 
ported a weaker dependence with the strongest by Hur et 
al.l^Sl being D ~ N^^-^. This is most probably due to the 
use of effectively shorter chains in the previous studies. 

Another surprising result regarding the diffusion of the 
rings is that the stress relaxation follows an approximate 
t~^^'^ scaling law, in fair agreement with the pre- 
diction of Kapnistos et al.'^ In fact, the TV = 800 ring 
system shows this behavior for at least four decades in 
time. While it is tempting to appeal to the Rouse model 
since it predicts G{t) ~ this certainly cannot ap- 

ply. This fast decay, compared to what we observe from 
the mean-square displacements of the beads, is in accord 
with the observation that Ccc(t) decays with a relaxation 
time Treiax ^ ]\[2.2±o.i _ That means stress can relax by 
internal reorganization without the need to diffuse a sig- 
nificant distance away from its original neighbor chains or 
to fully relax the overall shape of the chain. This however 
is only possible because the rings eventually are compact 
objects. As a consequence, the characteristic stress re- 
laxation time Troiax ~ A^^'^ implies that the very small 
viscosity 770 and the diffusion time, td ^ N^, x > are 
not directly related to each other. This line of arguments 
does not take into account the existence of, as we called 
them, quasi-entanglements, which relax only very slowly 
and might give rise to a plateau-like regime in the stress 
relaxation function. For a possible plateau, however, we 
expect a very small value of the order 0{kBT/N), simply 
because the number of neighbors that a given ring can 
entangle with becomes independent of N and settles at 
a value below 20, a typical number for a liquid. Thus, if 
existing, such a plateau will be very difficult to observe 
and beyond the accuracy of our data. 

An attempt to describe the amoeba-like motion of the 
individual ring among other nonconcatenated rings was 
presented in a recent article by Milner and Newhall^^. 
These authors argue that since the ring in the lattice of 
obstacles passes each lattice gate twice in the opposite 
directions and, therefore, represents a kind of tree, then 
each bond separates the whole system in two halves, with 
the sizes k and N — k, respectively. The smaller of the 
two numbers Milner and Newhall call "centrality" and 
provide convincing simulation data to the effect that ef- 
fective motion along the "reaction coordinate" of central- 



ity is sub-diffusive with Sk ^ t^^^. Clearly, the diffusion 
of centrality in this picture is also closely related to the 
stress relaxation, which means that this picture is ide- 
ologically very close to that proposed by Kapnistos et 
aP2l. Given the identical values of the exponents in the 
diffusion of centrality case and in our g^ (t) behavior for 
intermediate times and ring sizes, it is tempting to bridge 
the two statements. However, our conformations some- 
how suggest a rather extended regime of beads to qualify 
for the centrality argument, making the precise meaning 
of this concept somewhat unclear here. Future work will 
be required to understand whether they are really related 
or there is just a coincidence of numerical values for the 
two exponents. 



V. CONCLUSIONS 

The structural and dynamic properties of ring polymer 
melts have been investigated by molecular dynamics sim- 
ulation using a semiflexible bead-spring model in this and 
the previous paper. For the longest chain lengths, the ra- 
dius of gyration was found to approach a scaling regime 
of N^^^ in agreement with previous simulation results 
and theoretical studies. While this suggests collapsed 
conformations, the rings are, in fact, open or expanded 
and highly interpenetrated by their few 0{N^) neighbors. 
This is supported by a primitive path analysis which was 
also used to show that the estimated quasi-entanglement 
number of the rings is larger than that of the analogous 
linear chains. The rings being not as entangled as the 
linear chains and more compact in size may explain why 
they relax stress faster and have a viscosity which scales 
much more weakly than the N^''^ dependence of linear 
chains. Despite following different scaling laws with N 
for the radius of gyration, stress relaxation and viscos- 
ity, and undergoing different modes of motion while dif- 
fusing, for large N the two architectures approximately 
obey the same D ~ ^-2.3±o.i power law. Because of 
the compact nature of the ring conformations, resulting 
in a fairly deep correlation hole, obtaining diffusion con- 
stants from such simulations might be subject to rather 
strong finite size effects in the number of chains consid- 
ered. Such polymeric systems have some resemblance 
to a soft sphere liquid which is known to be subject to 
such effects. However, we are confident that our sys- 
tems with M = 200 ring polymers are sufficiently large 
and we emphasize that studies using significantly fewer 
chains should be avoided. 

Though we have presented a rather comprehensive 
analysis many questions remain open. Conformational 
properties need further attention, especially for questions 
with respect to biological issues. This concerns, for ex- 
ample, contact maps and the life time of the clusters 
seen in such mapa^. Also, the precise determination of 
(3 and 7 and the related scaling prediction /? + 7 = 2 
need further attention. The decoupling of translational 
motion and stress relaxation for rings, the latter leading 
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to the low viscosity, certainly is not fully understood and 
needs further investigation. The sub-diffusive regime in 
the MSD of the rings up to displacements that are sig- 
nificantly larger than poses special challenges for 
a proper theoretical description. To better compare our 
findings to experiments, a systematic study of the influ- 
ence of linear "contaminants" is needed since experiments 
never can exclude them completely. Such a study is in 
11 



preparation^ and will also include an extensive analy- 
sis of the dynamic scattering functions of both the rings 
and the linear polymers. Recent advances in synthesis, 
purification and characterization should allow a better 
comparison in the near future. While the main thrust 
of the two papers was on the properties of nonconcate- 
nated ring polymer melts, for comparison we also have 
presented extensive simulations of highly entangled melts 
of linear polymers. While in general the reptation pic- 
ture provides a good qualitative and semi-quantitative 
theoretical framework there are still open questions with 
respect to the slow crossover to the anticipated asymp- 
totic regime or possible deviations thereof. 
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